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We study in details the distribution of mass in galaxy cluster CL0024+1654 inferred using the 
method of strong gravitational lensing by Tyson et al. (1998). We show that a linear correlation 
exists between total, visible and dark matter distributions on log-log scale with consistent coeffi- 
cients. The shape and parameters of log-log-linear correlation are not affected significantly whether 
one uses projected or volume mass densities but is consistent with k = 2 — 5 visible/dark ratio. We 
also show and analyze in depth so called alignment properties of the above-mentioned profiles. We 
show that log-log-linear correlation and alignments can all be understood in terms of thermody- 
namic/hydrodynamic equilibrium with gravitational potential growing almost linearly in the region 
of interest. We then analyze the hypothesis of thermal equilibrium on the base of the existing data 
about CL0024 cluster. If the presence of log-log-linear correlation and alignments were interpreted 
thermodynamically, this would indicate the mass of the dark matter particle 2-5 times smaller than 
that of atomic hydrogen, thus giving range for the mass of dark matter particle between 200MeV 
£ . and lOOOMeV. 

<^- ' I. INTRODUCTION 

<n: 

In the past decade the field of astronomy and astrophysics has experienced a period of rapid development thanks 

to the operation of Hubble Space Telescope and a wide spread of new faster computers. This development brought 

(yT) ' a series of brilliant insights in the puzzles of the structure and history of the universe. At this time the growing 

\Q , number of precision measurements of galactic rotation curves and observations of gravitational lensing in galaxies 

l/"~) ' and galaxy clusters combined with extensive computer simulations shed new light on the problem of hidden, or dark, 

matter in the universe and its role in the universe evolution. Still, little is known about the microscopic properties 

and composition of the dark matter. Many possibilities have been put forward by various extensions of the Standard 

Model but yet existence of any of such particles is to be confirmed experimentally. 

Recently it had been pointed out that the distribution of dark and visible matter in spiral galaxies with various 
i-Q ' luminosities and galaxy clusters exhibit strikingly similar correlations with almost identical parameters [1]. Specifically, 
i*, it was found that in all of these systems the mass densities of visible and dark matter are linearly correlated on log- 
log scale. The proportionality coefficient in this correlation appears to be universal and equal to k ~ 3 — 5. This 
conclusion had been drawn from the analysis of dark and visible matter distribution in the galaxy cluster CL0024 [2] 
c/2 ■ and synthetic mass model for spiral galaxy rotation curves of Persic and Salucci [3] . 

In this paper we are going to further sharpen our focus on the properties of the mass distribution in galaxy cluster 
CL0024. We will investigate in details presence and properties of log-log-linear correlations, alignment properties and 
their thermodynamic significance and interpretation. 

We begin in the next section (Section II) with discussion of the properties of the projected mass density profiles 
S(_R) derived by Tyson et al. in 1998 and focus on the log-log- linear correlation between them. We then continue with 
analysis of the volume density profiles p{R) inferred from S(i?) via inverse Abel transformation and show that the 
log-log-linear correlation is present also between the volume density profiles. We comment on the significance of these 
correlations and their thermodynamic interpretation. Finally, we review the alignment properties of the mass profiles 
in the galaxy cluster CL0024, mentioned earlier in the literature, and show that they imply exponential behavior for 
the radial mass density £(i?) in the region of interest. In Section III we present analysis of the thermal state of the 
matter in the galaxy cluster based on the strong gravitational lensing study carried out by Tyson et al. Summary 
and conclusions follow in Section IV. 

II. MASS DISTRIBUTION IN THE GALAXY CLUSTER CL0024 

That the mass distribution in galaxies or galaxy clusters can be measured using gravitational lensing had been 
pointed out long time ago [4]. Gravitational lensing is one of the consequences of Einstein's theory of gravity (General 
Relativity) in which the light from a distant object is bent by a heavy galaxy or galaxy cluster to produce "fake" 
images of the original object. The amount of this distortion is very tiny and for an observable effect to be seen the 
objects' alignment should be quite perfect. Still, with billions of galaxies in the sky, a number of gravitational lenses 
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had been discovered in recent years providing excellent tests of the Theory of General Relativity. One of the most 
famous examples of gravitational lensing is the object known as Einstein's cross (or Q2237+030), in which the light 
from a distant quasar is disrupted by a low-redshift galaxy to produce four images in the form of a cross surrounding 
the galaxy's nucleus. 

Galaxy clusters are one of the most interesting objects for gravitational lensing studies because of their unsurpassed 
mass and large extent. In fact, precision measurement of gravitational lensing in a galaxy cluster is capable to produce 
a detailed map of gravitating mass distribution inside the cluster and thus provide valuable information about its 
structure. In the last decade many clusters exhibiting gravitational lensing have been found and their "gravitating" 
mass distributions have been analyzed [2,5-9]. It was generally shown that a large portion of mass in galaxy clusters 
is not associated with the luminous galaxies and form a smooth extended distribution. However, practically never 
such mass maps, obtained from lensing, have been compared with the distribution of the visible matter inside the 
corresponding clusters. The study of strong gravitational lensing in galaxy cluster CL0024 by Tyson et al. in 1998 
[2] is distinguished in that the detailed maps both for total and visible mass in the cluster were constructed and 
presented. 

The galaxy cluster CL0024+1654 is a remarkable instance of strong gravitational lensing in which multiple images 
of a background galaxy with distinctive spectrum are formed. In 1998, an analysis of Hubble telescope images of this 
cluster was carried out by Tyson et al. and the mass profile of the cluster was obtained from strong gravitational 
lensing. Tyson et al. found that the vast majority of the mass in CL0024 is not associated with the galaxies and forms 
a smooth elliptical distribution, slightly shallower than isothermal sphere, with a soft core of r core = 35 ± 3h~ 1 kpc, 
where h is the normalized Hubble constant. No evidence of in- falling massive clumps had been found for the dark 
component. The projected dark matter density profile was fit well by a power-law model 

*(y) = (1 + y2)2 -v W 

where y = r/r core , K = 7900 ± IOOIiMqpc^ 2 , r core = 35 ± 3h~ 1 kpc and r\ — 0.57 ± 0.02. The primary conclusion 
was that the cusped mass profile for the dark component, suggested by many-body simulations within the Cold Dark 
Matter model [10], was inconsistent with the observed results. Along with the total mass distribution Tyson et al. 
also presented the radial distribution of the visible matter density and of the visible light density. 

It was recently noticed that the mass profiles obtained in Ref. [2] possess interesting correlation properties [1]. 
Specifically, if the mass profile of the dark matter and that of the visible matter are plotted one vs. the other on 
log-log scale, the linear correlation between the two becomes apparent which implies 

log£„ w K^logSrf. (2) 

Remarkably, such behavior of mass profiles should be expected on thermodynamical grounds. For example, for 
isothermal distribution of self-gravitating gas it is known that the mass profiles satisfy 

logp^-e t , 
logp d -e t , 

where [ii is the molar mass for the corresponding component and $(r) is the gravitational potential at position 
r [11]. Then, independent from the details of the gas radial distribution, \ogp v ~ logp^. With somewhat more 
intricate calculation, the same conclusion can be drawn for a self-gravitating gas just in hydrodynamic and not in 
full thermal equilibrium. Consider a multi-component gas cloud in hydrodynamic equilibrium, then for the density of 
each component we can write 

^+*(r)^=0, (4) 

dr dr 

where Pi(r) — pj(r)Tj(r)//ij is the partial pressure for component i. After a simple manipulation we obtain 

T t {r) d log Pi (r) | 1 dT t {r) | rf$(r) 
Hi dr /i i dr dr 

and 



(5) 



\<KI> l (r) = C-i li I [dloglKr) + ^]. (6) 



This tells us that logp v ~ \ogp d ~ J(d\ogT + yrf$) in case the components have temperatures which are locally 
similar or equal. Thus, the above mentioned log-log-linear correlations may contain important information about the 
microscopic properties of the dark matter. In this paper we would like to investigate in more details the properties 
of the mass profiles in the galaxy cluster CL0024. Since in cosmology huge distances and small heat transfer rates 
seem to interfere with the processes of thermal equilibration, thus reducing the likelihood of thermal equilibrium, such 
study would be most certainly beneficial. 

In Rcf. [2] the projected radial density profiles for the total mass and the visible matter were originally presented 
on log-log scale and mass profile for the dark matter could be straightforwardly deduced from these. The total mass 
and the dark matter form smooth density distribution monotonously decreasing with distance. Similar behavior is 
observed for the visible matter except for a flat segment at distance of about 100h~ 1 kpc [see Fig.(l)]. As can be 
seen from Fig.(l), at this distance £„(i?) is approximately constant for about 35h~ 1 kpc while all the other profiles 
continue to fall. While we do not know the origin of such peculiar anomaly, we note that outside of this flat region 
the mass profile for visible matter has essentially the same behavior. In fact, if we cut out this flat segment and shift 
the tail of the visible matter distribution to make a continuous curve, we would observe that the obtained profile is 
almost a perfect straight line on logS — R scale [data shown with triangles in Fig.(l)]. Upon inspection, one notices 
that, actually, on logE — R scale all profiles fall linearly with the distance, which implies behavior Sj(i?) ~ e~ aiR . 
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The projected mass density profiles from Tyson et al. plotted on different scales. Up to down are total, dark 
matter densities, respectively; triangles correspond to distribution of visible matter corrected for the anomaly at 



As was already mentioned, it can be noticed that mass distributions of dark and visible matter in CL0024 exhibit 
linear correlation on log-log scale. Similar property can be observed also for other pairs of mass profiles. For example, 
if log of the total mass is plotted vs. log of the dark mass, the linear correlation becomes very prominent. The same 
is true for the total and the visible matter mass densities. In this latter case the anomalous flat segment becomes 
more prominent; however, if the visible matter mass profile is corrected for this region, the points on log-log scale 
form almost a perfect line [data shown with triangles in the diagram in Fig. (2)]. The relations between different mass 
profiles can be further analyzed by looking at different segments of the diagrams in Fig. (2) and Fig. (3). We found that 
the correlation coefficients tend to decrease somewhat if one moves away from the center of the cluster. We obtained 
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We also note that all three coefficients are consistent with each other and imply the concentration by mass of the 
dark matter about 80%. For example, from thermodynamic interpretation of Kij « (fii/fj,j)(Tj/Ti) we would expect 
KvtKtd ~ K vd and this is indeed the case. Furthermore, we can write for the total mass density 



M(r) = m d (r) + m v (r) rj a d e-^* (r)/T + a v e~»'*W r , 



(8) 



where md( v ) is density of dark (visible) matter and Hd(v) is t ne molar mass of dark (visible) matter. Then, for variation 
8M(r) we obtain 



6M(r) w m d (r)(-fi d 5$/T) + m v (r){-p v 8<$>/T) 



(9) 
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Utot 



m d (r) m v (r) 



M(r 



"McT 



M(r 



(10) 



is effective total molar mass. By dividing Eq.(lO) with fi v or ^ and taking into account that, for example, (Mot/Hd = 
Ktdi onc can immediately verify that all of the above-mentioned coefficients n satisfy Eq.(lO) with the dark matter 
concentration by mass nid/M ~ 0.8, consistent with direct observation from the data by Tyson et al. [see Fig. (4)] 
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FIG. 2. Linear correlations between mass profiles in galaxy cluster CL0024. Shown are total and dark matter (left panel), 
total and visible matter (center) and total and visible matter corrected for anomalous region (left panel, data shown with 
triangles) . 
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FIG. 3. Linear correlation between log of dark and log of visible matter mass density. Emphasized in the left panel is the 
anomalous behavior of the visible matter mass profile between 80h~ 1 kpc and 110/i -1 kpc. Also shown is the log- log- linear 
correlation between dark and visible matter corrected for anomalous region (data shown with triangles in the right panel). 
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FIG. 4. Concentration by mass of the dark (solid line) and the visible (dashed line) matter in CL0024 as function of distance 
R. Also shown is the ratio of these concentrations (thick dashed line) . 



The existence of log-log-linear correlation in CL0024 becomes even more striking once one realizes that similar 
relation exists in spiral galaxies. In spiral galaxies the presence of dark matter reveals itself through the anomalous 
behavior of the Rotation Curves (RC) which for large distances do not fall as 1/y/R, according to the Kepler's law, 
but remain roughly flat. Such behavior had been analyzed in the literature and fit by synthetic model of Persic and 
Salucci [3] in which a number of RC where categorized by galaxy luminosity and reproduced with superior accuracy 
with a simple analytic model (Universal RC) . In this model the distribution of visible matter in the galaxy is described 
by exponential thin disk with surface density 

/ t ,(r)^ e - 3 - 2r/ - R ^, (11) 

which is known to fit well the surface brightness of galaxy disks, and dark matter is represented with a spherical 



flat-core distribution. We noticed that the URC for all luminosities can also be successfully described by a spherical 
dark hallo with density 

p d (r)~e- a ' r / R °", (12) 

in which case log/„(r) and log pd{r) are linearly related with logi„(i?) w n vd log pd{R)- Surprisingly, the coefficient 
of such relation varies very insignificantly with galaxy luminosity and corresponds to n s J d tra ~ 2.6 — 3.8. This is 
consistent with k vc [ obtained from the analysis of log-log-lincar correlation between the mass profiles in the galaxy 
cluster CL0024. Considering that such log-log-linear correlation is, in fact, expected from thermo- or hydrodynamic 
equilibrium, the above observation is most indicative of the presence of a sort of equilibrium between the dark and 
the visible component in spiral galaxies and galaxy clusters. 

Let us now concentrate on the volume densities of visible and dark matter in CL0024. Note that the log-log-linear 
correlation observed above was obtained for the projected mass densities £j(i2), which are the quantities actually 
measured in the gravitational lcnsing. At the same time Pi(R), appearing in Eq.(3), is the volume mass density related 
to T,i(R) via relation 

E(r,0) = dzp(r,6,z). (13) 

The projected mass density in general does not uniquely specify the volume density. Only with additional constraints, 
e.g. the spherical symmetry, Eq.(13) can be inverted; in spherically symmetric case this procedure is known as inverse 
Abel transformation [12]. As is shown elsewhere, if mass distribution p(R) decreases fast enough, the projected mass 
density can be approximately represented as 

E(r) « h(r)p(r), (14) 

where h(r) is a slow varying (h(r) ~r") effective depth and p(r) is the density maximum along the line of integration 
in Eq.(14) [13]. The correlation between log £i(i?) then implies 

log Si (r) « on log r + log p t (r) ~ ^ 

Kij log T,j (r ) w Kij (oij log r + log pj(r)) 

and, if log p(r) varies quicker than logr, we get linear correspondence 

log ^ (r) n Kij log pj(r). (16) 

Since we have observed before that the mass profiles Si (R) in the galaxy cluster CL0024 decrease exponentially with 
distance, we do expect that linear correlation will be present between p's as well. We can estimate the error introduced 
in our treatment because of the use of projected mass density S(i?) in place of the volume density p(R) by inverting 
directly the experimental profiles £(i?). 

Unfortunately, the inverse Abel transformation requires knowledge of £(i?) for arbitrary large R while ex- 
perimentally we only know a limited segment R < R max ~ 120h kpc for visible and dark components and 
R < Rmax ~ 220h~ 1 kpc for total mass profile. Certain strategy must be used to extrapolate the experimental 
profiles to R > R ma x- In this study we adopt three main strategies: cutting the integration in the inverse Abel trans- 
formation to R < Rmax, thus taking E(i? > R m ax) = 0; fitting profiles with a sum of exponentials S(i?) = ^2 a,ie~ R ' Ri 
or fitting profiles with a sum of power-law functions S(i?) = y^qj(l + R/Ri) hi . As should be expected, in every ap- 
proach obtained radial profiles are similar for R small; however, at larger distances a significant model-dependence 
becomes increasingly apparent. In fact, obtained mass profiles are only model- independent for the central region of 
the galaxy cluster out to about R m ax/2 [see Fig. (5)]. At the same time the central part of these profiles is most 
sensitive to the noise in S(i?), so that we inevitably obtain large errors and noise in p(R) practically for all values of 
R [see Fig. (6)]. 




FIG. 5. Model dependence in volume density profiles obtained by inverse Abel transformation. Solid line is extrapolation 
E(7? > Rmax) — 0, doted line is exponential fit and dashed line is power-law fit. 
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FIG. 6. Radial volume density profiles obtained by inverse Abel transformation for total mass, dark matter and visible 
matter (left, center, right respectively). Points correspond to S(_R > R ma ,x) — extrapolation, dashed line corresponds to 
exponential fit and thick dotted line corresponds to power-law fit. 
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FIG. 7. Inverted radial mass densities compared on log-log scale. As before, triangles denote the visible matter distribution 
corrected for the anomalous region. 



While data-points after inverse Abel transformation are more noisy, the linear correlation still can be seen on the 
log- log scale between any two mass profiles [see Fig. (7)]. The parameters of this correlation are somewhat dependent 
on the way the original S-distributions were extrapolated, but such variation is well within acceptable bounds. Taking 
into account this variation, for the parameters of log-log-lincar correlation between volume density profiles we obtain 
K v t ~ 1.2 — 2.1; Ktd ~ 1.0 — 2.0 and n v d ~ 2.1 — 3.4. As we can see, the correlation coefficients are close and generally 
somewhat less from those derived from the projected mass densities; the average error introduced using the projected 
mass density S(-R) in place of volume mass density p(R) is of the order of 30%. 

To understand the origin and magnitude of this error, let's consider effective depth defined by S(i?) = h(R)p(R). 
Then 



logZ(R) = log h(R) + log piR), 



(17) 



and if p{R) falls off rapidly with R, one expects that log h(R) will vary relatively slowly with distance and thus 
the linear correlation between log£(i?) would imply a linear correlation between logp(i?). Nonetheless, variation of 
log h(R) would affect the estimate for k and introduce an error of the order of A log h(R)/A log p(R) . For example, for 



the profile of dark matter in the region of interest A log h s=s 0.5 and A log p ; 
by this in our estimate of n v d is of the order of A log h/A log p ss 30%. 



1.5. As is expected, the bias introduced 




Log Y. 



Log R 




Log R 



FIG. 8. Alignment of the projected mass profiles in CL0024 by log-scaling (left) and distance-scaling (right). 
Let us now turn our attention to the alignment properties of the mass profiles in CL0024. In Ref. [1] it was noticed 



that all mass profiles can be aligned very well by either log-scaling [logS,(r 



logSj(r)] or distance-scaling 



[log Si (r) <~ log Sj (re~ Zii )] , as can be seen from Fig. (8). For the parameters of such alignment we find that the dark 
matter profile can be aligned with the total mass profile by rescaling with Ktd ~ 1.25 — 1.45, the mass profile for 
visible matter can be aligned with that for total mass with n vt <~ 2.0 — 2.5 and the visible mass density profile can be 
aligned with that for dark matter with k vc i ~ 2.5 — 4. These values are similar to those obtained from log-log-lincar 
correlation. In fact, it becomes obvious that log-scaling alignment is directly related to the existence of log-log-lincar 
correlation and vice-versa. 

Similarly, all profiles can be aligned by distance-rescaling S(r) — > S(r • e~ z ). We found that the best fit parameters 

are z v d ~ —0.15 0.25 for visible-to-dark alignment, z vt « 0.45 — 0.5 for visible-to-total alignment and Zdt ~ 0.6 — 0.9 

for dark-to-total alignment. The implication of these alignment properties is less obvious; however, it can be shown 
that together with log-scaling they imply 



log£*(r) 



,107(log(r)-*i) + hi 
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(18) 



From the parameters k and z listed above we find 7 w 0.8 — 1.0. This conclusion is consistent with our earlier 
observation that the projected mass density in galaxy cluster CL0024 fall exponentially with distance; 



S i (r) = B i e-°( r / ri ) 7 ,7«l- 



(19) 



If this result to be interpreted thermodynamically, log p(R) ~ $(-R), that would imply that the gravitational potential 
is rising in the region of interest almost linearly. Let us note that the same alignment properties are observed for the 
volume mass densities as well. All three volume density profiles can be aligned very well with the above-mentioned 
scaling transformations [see Fig. (9)]. The parameters of these alignments are Ktd ~ 1-5 — 1.8, K vt s=s 1.7 — 1.9, 
K v d « 2.8 — 3.4 and z v d ~ —0.35 0.5, z vt s=s 0.4 — 0.7 and Zdt ~ 0.6 — 0.8, similar to those obtained from T,(R). 
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FIG. 9. Alignment of the volume mass density profiles in CL0024 by log-scaling (left) and distance- scaling (right). 



III. MASS PROFILES AND THE THERMAL STATE OF MATTER 

As was mentioned above, the log-log-linear correlation between mass profiles of the galaxy cluster CL0024 is a 
very remarkable property, especially given that a similar observation holds for a variety of spiral galaxies and that 
such correlation is actually expected in case of thermodynamic or hydrodynamic equilibrium. The most tempting 
interpretation of such correlation is through the thermo- (hydro-) dynamic equilibrium, in which case the k vc [ rj 3.0 
should be related to the ratio of the molar masses for visible and dark matter [1]. This assertion, however, is difficult 
to understand theoretically from what is known or assumed today about the properties of the dark matter [13]. While 
there exist no direct data indicative of the thermal state of the dark matter, the commonly accepted notion is that 
the dark matter is extremely weakly interacting and, thus, may well be thermally isolated from the visible matter. 
Here we would like to see if one can obtain any information about the thermal state of the matter in galaxy cluster 
CL0024 from the actual mass profiles provided by Tyson et al.. 

According to Eq.(3) and Eq.(4), the temperature of matter may be estimated from p(R) using 



(j, 



log p(R)- log p(0) = ^($(0)-$(R)) 



(20) 



or from solving the differential Eq.(5) with respect to T(R) with p(R) known. In the latter approach one needs to 
overcome the difficulty that the result of integration in Eq.(5) depends on unknown integration constant, for example 
T{R max ). This problem may be partially avoided by integrating Eq.(5) from large distance toward smaller r to reduce 
dependence on the unknown boundary condition T{R max ) or by using T(R max ) close to that obtained from Eq.(20). 
In either case it can be observed that in the region of interest the dependence on boundary condition T{R max ) is 
weak. 
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FIG. 10. Temperature profile extracted from £(i?). In the left panel shown is the temperature extracted from Eq.(20) and in 
the right panel shown is the temperature extracted from Eq.(5). Triangles (higher line) correspond to dark matter temperature 
profile, boxes (lower line) correspond to visible matter and diamonds correspond to effective temperature extracted from the 
total mass distribution. 



If the above idea is applied to the galaxy cluster CL0024 with the assumption p(R) ~ £(i?), as above, the temper- 
ature extracted in this way indicates that the dark matter has smooth temperature profile cooling outward, while the 
visible matter is mostly isothermal with temperature somewhat decreasing in the center of the galaxy cluster. The 
estimate for the temperature of the visible matter T v s=s 5 • 10 8 K and its decrease toward the center of the cluster is 
consistent with our understanding of X-ray clusters and the effects of radiative cooling in their centers [14]. Results 
obtained in this way, however, are only qualitatively reliable. 

Since logp(i?) « — log h(R) + log£(i?), in the case when the temperature profile T(R) is almost isothermal even 
small correction log h(R) may be significant. Generally, we expect the "effective" depth log h(r) to increase the 
temperature of the distribution £(i?) in the center of the cluster relative to the temperature associated with p(R). 
Obviously, that may have a dramatic effect on our results making, for example, the temperature profile for the dark 
matter isothermal or breaking isothermality of the visible matter. 

Thus, we do need to consider the above approach applied directly to the volume density profiles p(R). As we have 
noticed before, the data points after this transformation are significantly more noisy. This problem is even more 
serious here than in the case of log-log-lincar correlations since the resultant T(i?)-profile is particularly sensitive to 
the noise in p(R). In fact, the volume density profile obtained with S(i? > R ma x) — is so noisy that hardly any 
conclusion can be obtained from it at all. As can be seen in Fig. (11), the only conclusion we can draw is that the 
temperature of the visible matter is definitely decreasing toward the center of the cluster and that the "effective" 
temperature obtained from the total distribution is decreasing in the outer regions of the galaxy cluster. 
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FIG. 11. Temperature profiles from volume densities deprojected with T,(R > R max ) = 0. Triangles denote dark matter 
temperature, boxes correspond to visible matter and diamonds correspond to effective temperature of the total mass profile. 

Better results can be obtained when the mass profiles S(i?) arc smoothencd and extrapolated with a linear combi- 
nation of exponential or power-law functions. For both fits we obtain a similar result that the visible matter is cooling 
significantly in the center of the cluster, consistent with the strong radiative cooling in X-ray clusters mentioned 
above. At the same time for the dark matter we obtain the distribution close to "isothermal" which typically cools 
very slightly in the center. It becomes obvious, therefore, that the local thermal equilibrium between dark and visible 
matter is significantly broken, at least, in the center of the galaxy cluster. Unfortunately, while these qualitative 
conclusions are similar for both exponential and power-law extrapolations, the quantitative details differ dramati- 
cally between the two, especially for the outer regions of the cluster. The integration of Eq.(5) using exponential 
extrapolation produces temperature distributions for the dark and the visible components with almost constant ratio 
K v d ~ 1.8. This ratio gets larger toward the center of the galaxy cluster as the visible matter cools rapidly and 
the dark matter fails to follow the suit. The result of analysis for this case, thus, supports the hypothesis of overall 
thermal equilibrium. It implies that, although the ratios T/fi for visible and dark matter vary with distance, the ratio 
K v d = (Td/T v )(/j, v / fid) remains constant for a wide range of values of R. The thermal equilibrium is only broken in 
the central part of the cluster where the rapid radiative cooling of the visible component is essential. 

However, the conclusion obtained using the power-law extrapolation is exactly the opposite. In this approach we do 
find that at large distances the ratios T//J, for the visible and the dark matter is practically the same with k v< x rj 1. This 
would be consistent with our understanding of the cluster formation. According to our current understanding, the 
primary heat source in galaxy clusters is gravitational heating. Then, as originally cold gas collapses thereby forming 
galaxy cluster, the gas molecules acquire kinetic energy AEi s=a /iiA$, where A$ is the change in gravitational 
potential during the collapse. Later this kinetic energy is transfered to thermal energy of the matter inside the 
cluster, so that Tj s=s ^ A$ and thus Tj//ij s=s A$ is the same for all kinds of the particles in the gas cloud. We expect 
therefore the mass distributions for different components to be similar and k vc i ~ 1 , consistent with the result obtained 
with power-law extrapolation, k is not equal to 1 only in the center of the cluster where the visible distribution cools 
rapidly while the dark matter distribution remains almost isothermal. The temperature of the visible matter in both 
approaches comes out T v s=a 10 8 K. 
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FIG. 12. Temperature distributions for exponentially extrapolated E(_R). Shown are T(R) given by Eq.(20) (left panel), 
obtained from Eq.(5) (central panel) and extrapolation to large distances (right panel). Higher curve corresponds to dark 
matter, while lower curve corresponds to visible matter and the middle curve corresponds to "effective" temperature obtained 
for total mass profile. 
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FIG. 13. Temperature distributions for power-law extrapolated T,(R). Shown are T(R) given by Eq.(20) (left panel), obtained 
from Eq.(5) (central panel) and extrapolation to large distances (right panel). Higher curve corresponds to dark matter, while 
lower curve corresponds to visible matter and the middle curve corresponds to "effective" temperature obtained for total mass 
profile. 

We must conclude, therefore, that the data on the distribution of mass in CL0024 are not sufficient neither to 
support nor to reject the hypothesis of thermal equilibrium. Only general conclusion can be drawn from this analysis. 
In particular, the temperature extracted from the mass profile of the visible matter is of the order of 5 • 10 7 — 1.5 • 10 & K 
and is consistent with X-ray data on galaxy clusters. The temperature of the visible component drops significantly 
in the center of the cluster and also falls slowly at large distances. The dark matter distribution appears to be close 
to isothermal with indications of slight cooling in the center. We believe that such behavior may be indicative of 
nonvanishing thermal connection between the dark and the visible components in the cluster. Still, the data are too 
uncertain to make this assertion any more definite. 

The assumption of the local thermal equilibrium is most certainly violated, at least in the central region, where the 
effect of radiative cooling of the visible component is the strongest. In the other parts of the galaxy cluster, however, 
we can neither confirm nor reject the hypothesis of thermal equilibrium as conclusion depends strongly on the way to 
interpolate the available data to large values of R. 

IV. CONCLUSION 

In this work we extend the analysis of the visible matter, dark matter and total projected mass profiles derived for 
the galaxy cluster CL0024 from strong gravitational lensing by Tyson et al.. We observe that the linear correlation 
exists between each pair of mass profiles on log-log scale. This linear correlation is preserved whether one uses the 
projected mass profiles [£(./?)] or the volume density profiles [p(i?)]. The difference between the coefficients of log- 
log-linear correlations obtained for E(i?) and p(R) is small and is of the order of 30%. We obtained following values 
for the parameters of log-log-linear correlation between dark and visible mass distribution profiles 



K v d ~ 2.9 — 4.4 (from projected mass profiles £), 
k vc [ ~ 2.1 — 3.4 (from volume density profiles p), 
k vc i ~ 2.2 — 4.0 (from alignment properties for E), 
k vc i ~ 2.8 — 3.3 (from alignment properties for p). 



(21) 



The correlation coefficients for other pairs of mass profiles were also listed in this paper. We also analyzed in details 
the alignment properties of the mass profiles in the galaxy cluster. We found that the log-scaling and the distance- 
scaling, mentioned earlier in the literature, imply (and are consistent with) exponential behavior of the mass profiles 
in the region of interest 



-a{r/rO)~ l 



7 



0.8-1.0. 



(22) 



These properties of the mass profiles are striking, especially since similar correlations are observed in the spiral 
galaxies with n v d ~ 2.6 — 3.8. Given huge differences between the systems where such correlations are observed, it 
appears that this cannot have an accidental character but is related to actual properties of dark and visible matter. 
Interestingly, such behavior of the mass profiles in self-gravitating system is expected in the case of thermo- or 
hydrodynamic equilibrium. Whenever the thermal states of the components in such system are similar, one expects 



\ogp v ~log/9d- 



(23) 



It is therefore tempting to interpret the observed correlations as a sign of equilibrium between the dark and the visible 
components in the galaxy cluster in which case the ratio K v d corresponds to the molar mass ratio of visible and dark 
matter 



10 



^w 2.1 -4.4. (24) 

This suggests the mass of the dark matter particle between fid ~ 200 — lOOOMeV. To our knowledge, there are 
no candidates within this mass range in the current extensions of Standard Model. Massive neutrinos and axions 
have experimental mass limits < 25MeV for the heaviest r-neutrino [15] which is well below our range. The most 
favorable candidates for non-baryonic dark matter, such as Weakly Interacting Massive Particle (WIMP) with the 
mass anywhere between lOGcV and ITeV and the SUSY lightest particle, e.g. neutralino, with the mass above 
30GeV [16], are far above the range obtained in our study. In fact, the closeness of fid to QCD-energy scale Aqcd 
may indicate that dark matter is ultimately related to QCD phenomena, such as quark and gluon condensations. 
Whether this is true and what kind of connection may exist between the dark matter and QCD is an interesting topic 
for further discussion. 

The concept of thermo- (hydro-) dynamic equilibrium between dark and visible matter encounters significant the- 
oretical difficulties in trying to understand the possible source of equilibrium between, supposedly, extremely weakly 
interacting dark and visible components. For that reason we further analyzed what kind of information about the 
thermal properties of the matter in the galaxy cluster CL0024 can be extracted from available measurement. In our 
analysis we found that for visible matter the average temperature is of the order of T v w 10 8 iG consistent with what 
we know about the temperature of the intergalactic gas in X-ray clusters. We also found that the temperature of 
the visible matter drops significantly in the center of the cluster, also consistent with what we know about radiative 
cooling in X-ray clusters. For the dark matter we obtained almost isothermal temperature profile with signs of slight 
cooling in the center of the cluster. We believe that such behavior may indicate existence of thermal connection be- 
tween dark and visible matter. At the same time we found that the available data are insufficient to either confirm or 
reject hypothesis of thermal equilibrium as the conclusion depends strongly on the way the data are smoothened and 
extrapolated at large distances. While in the central part of the cluster the equilibrium is most certainly broken by 
rapid radiative cooling of the visible component, in the other parts of the galaxy cluster one can find signs consistent 
either with full thermal equilibrium or full thermal isolation between the dark and the visible components - depending 
on the way the projected mass density profiles are extrapolated. 

We must conclude therefore that no reliable quantitative conclusion about the thermal state of matter in the galaxy 
cluster CL0024 can be reached on the basis of known mass density profiles. 

The observed similarities in the relative behavior of dark and visible matter in spiral galaxies and galaxy cluster 
CL0024 are very puzzling and deserve further investigations. While the precision of individual mass distribution 
measurements in galaxy clusters may be insufficient to obtain certain conclusions about the origin of such behavior, 
the synthetic analysis of few such studies may be helpful. One needs to see if the log-log-linear correlation can be 
found in other galaxy clusters and if its coefficients have similar values to those we have observed here. 
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